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Abstract 

Arbitrary amplitude ion-acoustic solitary waves (lASWs) are studied using Sagdeev-Potential 
approach in electron-positron-ion plasma with ultra-relativistic or non-relativistic degenerate 
electrons and positrons and the matching criteria of existence of such solitary waves are numerically 
investigated. It has been shown that the relativistic degeneracy of electrons and positrons has 
significant effects on the amplitude and the Mach-number range of lASWs. Also it is remarked that 
only compressive lASWs can propagate in both non-relativistic and ultra-relativistic degenerate 
plasmas. 
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I. INTRODUCTION 



Of the nonlinear excitations, ion-acoustic solitary waves (lASWs) are of the most im- 
portant and well-understood characteristics of plasma environments. Theoretical studies of 
main properties of these solitary structures date back to 1961 using Sagdeev pseud-potentials 
method [l|. Another method which is widely used to investigate the collective wave phe- 
nomenon in plasma is the so-called multi-scales perturbation method {2 8|. However, the 
latter method, which is based on approximation, is used only for the small-amplitude treat- 
ment of plasma in a state away from thermodynamic equilibrium. Therefore to obtain 
a good agreement with experiments, in this method, one needs to take higher-orders in 
perturbation amplitudes. In recent years there have been many investigations on solitary 
lASWs as well as solitary electrostatic waves (ESWs) in diverse plasma environments using 
Sagdeev pseudo-potential approach j9-13|. The small amplitude propagation and interac- 
tion of lASWs with relativistic degeneracy pressure effects have been recently considered in 
Ref. (1 

Among different kinds, pair-plasmas have attracted special attention in recent years, a 
special cases of which can be electron-positron (EP) and electron-positron- ion (EPI) 15l-ll9| 



plasmas. Electron-positron- ion plasma exists in places such as active galactic nuclei |20| . 
pulsar magnetospheres j2l| and in many dense astronomical environments, namely, neutron 
stars and white dwarfs [22] and may play a key role in understanding the beginning and 
evolution of our entire universe 23|. This kind of plasma may also be practically produced 



in laboratories 
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271]. More specifically when positrons, due to their significant lifetimes. 



are used to probe particle transport in Tokamaks, two component electron-ion (e-i) plasma 
behaves as three component (e-p-i) plasma 28|. Furthermore, the wave properties such as 
stabilities of a two component electron-ion (EI) plasma solitary excitations may be radically 
altered by inclusion of low amounts of positrons. 

Owing to their wide applicabilities in micro- and nano-electronic miniaturization, dense- 
plasmas is becomingone of the interesting fields of theoretical as well as experimental fields 



of plasma research 
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35| . Dense plasma or the so-called quantum plasma are characterized 



by high densities and low temperatures, however, a dense plasma may be realized in such 
hot places as planet interiors and white dwarfs [361] . More recently, quantum hydrodynamics 
model has been applied to study the electron-hole dynamics in semiconductors [37|, 
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The quantum effects in collective behavior of a plasma system becomes important when 
the inter-particle distances are comparable or less than the de Broglie thermal wavelength 
= h/(2nm(.kBTy^^ or equivalently when the thermal energies of plasma species are less 
than Fermi-energies [39]. In such cases the plasma becomes degenerate, in which the plasma 
ingredients are under effective influence of Pauli exclusion principle and classical statistical 
assumptions break down. Quantum effects also play important role in the nonlinear processes 
of white-dwarfs For instance, for a cold neutron star the densities can be as high as 
10^^^gm/cm~^ in the core, which is several times the density of an atomic nuclei. In extreme 
conditions such as the middle of a supernova or the core of a massive white dwarf the 
densities can be even catastrophically higher. At these very high densities the electrons and 
positrons may become ultra-relativistic giving rise to the collapse of star under its giant 



The present study is devoted to investigation of propagation of lASWs in an unmagnetized 
EPI plasma using Sagdeev pseudo-potential method in such extreme condition, taking into 
account the relativistic degeneracy effects for electrons and positrons. The basic normalized 
plasma equations are introduced in section [Tll Nonlinear arbitrary-amplitude solution is 
derived in section IIIII Section IIVI devotes to short argument about small amplitude lASWs. 
Numerical analysis and discussion is given in section |V] and final remarks are presented in 
section |Vll 

II. BASIC PLASMA STATE EQUATIONS 

Consider a dense plasma consisting of electrons, positrons and positive-ions. Also, sup- 
pose that the electrons and positrons follow the zero-temperature Fermi-gas statistics, while, 
ions behave as classical fluid. In such plasma electrons and positrons may be considered 
collision-less due to Fermi-blocking process caused by pauli exclusion principle. Therefore, 
the semi-classical description of nonlinear dynamics and interaction of waves in such plasma 
can be studied in the framework of conventional hydrodynamics model. The basic normal- 
ized equations describing plasma dynamical state may be written as 



gravitational force [41, 




where, electrons and positrons are of Thomas-Fermi type 

Tie = (1 + ^)hnp = a{l - aF^)^,a = ^,aF = (2) 

^eO -L Fp 

In obtaining the normahzed set of equations following scalings are used 

Vpe- , t - (0) - ^"^ksTpe , . 

X — )■ X, t —J- — , n — )■ nn^ % u — )■ uvpe, ^ — ^ ^ , (oj 

OJpi ^pi 6 



where, Upi = y e^rSe'^ / SQirni and vpe = ^/'ikBTpe/^i are characteristic plasma-frequency 
and electron Fermi-speed, respectively, and ni°'* denotes the equilibrium electrons density 
(ni^"* = ^p^pe with ppe being the electron linear Fermi-momentum). In a fully degenerate 



Fermi gas one may write the electron degeneracy pressure in the following general form 41 1 
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3/i3 



(2r2 - 3) Vl + r2 + 3sinh-Vl , (4) 



where, h and c are Plank constant and light-speed, respectively, and r = ppe/rrieC is the 
normalized relativity parameter. The electron number density can then be defined in terms 
of the relativity parameter {n^ = ^"^i^ r^)- It is noted that in the limits of very small and 
very large values of the relativity parameter we obtain 

^(f ) ^hcni (r — )► oo) 

Therefore, in a three-dimensional non-relativistic zero-temperature Fermi-gas for degenerate 
electrons and positrons from standard definitions we obtain Epj = (j = e,p) or 

Epj oc n^/o , which follows that ap = a"^/^. On the other hand, three-dimensional ultra- 
relativistic Fermi-gas, we have Epj = chkpj or ap = a^^^^. It is noted that in our model the 
inertial ions are always non-relativistic, hence, the non-relativistic hydrodynamics equation 
has been used in Eqs. ([T]). Therefore, the Poisson equation reads as 

= {l + ip)2 -a{l-apip)^ -rii. (6) 

At the equilibrium situation the overall neutrality condition gives rise to the following rela- 
tion 

/3 = l-a, /3 = !^. (7) 
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III. ONE-DIMENSIONAL ARBITRARY- AMPLITUDE ANALYSIS 



In this section we derive an appropriate Sagdeev pseudo-potential describing the dynamics 
of arbitrary-amphtude lAWs in plasma containing classical heavy positive ions and inertial- 
less relativistic or non-relativistic Thomas-Fermi electrons and positrons, obeying the three- 
dimensional distributions in Eqs. (|2]). Using one-dimensional version of Eqs. ([1]), in a 
reduced coordinate rj = x — Mt (M being Mach number which is a measure of soliton speed 
relative to ion-sound speed), from continuity and momentum equations we obtain 

1 



a 



2ip 



(8) 



where, we have used the fact that — )■ 0, -Uj — and Ui ^ f3 at rj ^ ±oo. Now, substituting 
Eq. ([8]) and Eq. (|2]) in Poisson's equation in Eq. ([1]), multiplying by ^ and integrating 

I — 1*. n fnr m — 1». -l-rsr^ \TiTP> Ai^r\\ri^ 

drj - 

2 



with boundary conditions {(f, ^} — )■ for 77 — ±00, we derive 



0, 



(9) 



where, the Sagdeev pseudo-potential V{ip) reads as 



2 r 



l-fl 



\5/2 



+ 



2a r 



5(7; 



1 - (1 - ap^ 



1 - 



M2 



(10) 



For the reality of V{(f) to be ensured we must have Lf < and 9? < ap^. The possibility 

of lAWs, therefore, require that 



dV{ip) 



dip 



0. 



tp=0 



dip^ 



<0, 



(11) 



ip=0 



it is further required that a exists such that V{ipm) = and for every ipm > <p > then 
V{ip) < 0. 

IV. SMALL AMPLITUDE THEORY 



Let us consider the small-amplitude limit in the above analysis. Expanding the potential 
V{ip) in (ITOl) near (p = 0, we obtain 



V" V" 



(12) 
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where, V^' = = 0) and V^" = V"\ip = 0) are computed from Eq. ([TOD as 



, _ 3 3ao> 1-a 



3 3aa| 3(l-a) 

^° " 4+ 4a + • ^^^^ 

Inserting into EqQ and integrating, we obtain (provided that Vq < ) a sohtary solution 
in the form (see ^]) 



This pulse profile is identical to the soliton solution of the Korteweg-de Vries (KdV) 
which is obtained by use of the reductive perturbation method, for example see 



equation, 



4^. It is 



important to notice that the soliton width L = 2/a/— Vq" and amplitude (po = —'^Vq/Vq" 
satisfy ^poL"^ = 12/Vq" = est., as known from the KdV theory. 

V. NUMERICAL ANALYSIS AND DISCUSSION 

As it was pointed out in section llllt the ion-acoustic solitary waves (lASWs) exist if the 
Sagdeev pseudo-potential satisfies the following conditions: (i) d'^V{{p)/dip'^\^^Q < 0, which 
reveals that the fixed point (f = is unstable; (ii) V{(fm) = 0, where, is the maximum 
value of ip; and (iii) V{ip) < when ipm > ip > 0. 

Noticing these conditions. Fig. (la) and Fig. (lb) show the areas in M-a plane, where: 
namely, lA solitary waves can exist for non-relativistic and ultra-relativistic electron-positron 
degeneracy, respectively. It is remarked that, the minimum values of Mach number, M, 
decreases as the fractional positron to electron number-density ratio, a, increases for ultra- 
relativistic case (see Fig. 1(b)). Nevertheless, the minimum values of M increases and 
reaches to a given maximum value then decreases as a increases (see Fig. 1(a)). On the 
other hand, the maximum value of M increases as a increases up to the value a ~ 0.37 
(a ~ 0.28) for non-relativistic (ultra-relativistic) case. The maximum value of M decreases 
as a increases in the range 1 > a > 0.37 (1 > a > 0.28) for non-relativistic (ultra-relativistic) 
case. 
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Comparing Fig. (la) and Fig. (lb) reveals that both supersonic and subsonic lASWs 
can propagate in ultra-relativistic case, whereas, only subsonic propagations can occur for 
non-relativistic case. Although both of supersonic and subsonic lASWs can propagate in 
ultra-relativistic degeneracy case, however, the former exist only for very small range of a, 
for approximately, 0.07 > a > 0.4. For a = 1 i.e., in the absence of ions, V{lp) does not 
depend on Mach number, M. This case corresponds to the solution of Poisson equation. 

In Fig. (Ic) and Fig. (Id), we have numerically analyzed the Sagdeev potential (Eq. fillip ) 
and investigated the effects of allowed values of M and a on the profile of the potential-well 
for both cases of non-relativistic and ultra-relativistic degeneracy. It is remarked that for 
fixed a (M) value, the increase of Much number M {a values) gives rise to an increase of 
both the potential depth and amplitude. The energy equation Eq. (9) has been numerically 
solved for some values of a and M. The potential pulse profiles has been depicted in Figs. 
(2a-2c). It is obvious that the potential profile becomes taller and narrower by increasing 
a and M for both types of non-relativistic and ultra-relativistic degeneracy, which is in 
agreement with above result. Also, we note that higher pulses are narrower, while shorter 
are wider, in agreement with the existing soliton phenomenology. Another important result 
is that the amplitude of pulse in non-relativistic case is higher than ultra-relativistic one for 
fixed values of a and M (Fig. 2(c)). Finally, we note that rarefactive lASWs do not exist 
for both of non-relativistic and ultra-relativistic degeneracy cases. 



VI. CONCLUSION 



The Sagdeev-Potential approach was used to investigate the propagation of ion-acoustic 
solitary waves (lASWs) in electron-positron-ion plasma with ultra-relativistic or non- 
relativistic degenerate electrons and positrons. The matching criteria of existence of such 
solitary wave were numerically investigated for both cases of ultra-relativistic or non- 
relativistic degenerate electrons and positrons. It is remarked that the characteristics of 
nonlinear lASW propagation differ in the mentioned cases. Both supersonic and subsonic 
lASWs can propagate in ultra-relativistic case, whereas, only subsonic propagations can 
occur for non-relativistic case. Only compressive lASWs can propagate in both of de- 
sired non-relativistic and ultra-relativistic degenerate plasmas. Also, it was concluded that, 
the differences tend to amplify by moving towards smaller values of fractional positron-to- 
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electron number-density a. In this work we consider cold ions and the effects of warm ions 
on lASWs in such plasmas can be investigated in future. 
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FIGURE CAPTIONS 



Figure- 1 

(Color online) The stability regions (dark) of arbitrary amplitude lASWs is shown in a — 
M plane for non-relativistic (Fig. 1(a)) and ultra-relativistic (Fig. 1(b)) electron/positron 
degeneracies. Figures 1(c) and 1(d) show the corresponding pseudo-potential dips for varied 
fractional positron to electron number-densities, a = 0.3 (blue), a = 0.4 (grey) and a = 0.5 
(red), and fixed Mach-number in non-relativistic (Fig. 1(c)) and ultra-relativistic (Fig. 1(d)) 
electron/positron degeneracies, respectively. 

Figure-2 

(Color online) (a) Ion acoustic solitary wave profiles obtained by numerical solutions of 
Eq. iQ corresponding to pseudo-potentials shown in Fig. 1(d) which correspond to the 
values a = 0.25 (black), a = 0.5 (blue) and a = 0.75 (red), (b) Profiles for different 
supersonic values of Mach-number M = 1 (black), M = 1.05 (blue) and M = 1.1 (red), (c) 
Comparison between non-relativistic (thin-line in black) and ultra-relativistic (thick-line in 
blue) electron/positron degeneracies for similar values of a and M. 
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g=0.25,0.5,0.75, M=0.8 (Ultra-relativistic) a=0.3, M=l, 1.05, 1.1 (Ultra-relativistic) 
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